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We calculate numerically the quasiparticle effective mass (m*) renormalization as a function 
of temperature and electron density in two- and three-dimensional electron systems with long- 
range Coulomb interaction. In two dimensions, the leading temperature correction is linear and 
positive with the slope being a universal density independent number in the high density limit. We 
predict an enhancement of the effective mass at low temperatures and a non-monotonic temperature 
dependence at higher temperatures {T/Tf ~ 0.1) with the peak shifting toward higher temperatures 
as density decreases. In three dimensions, we find that the effective mass temperature dependence 
is nonlinear and non-universal, and depends on the electron density in a complicated way. At very 
high densities, the leading correction is positive, while at lower densities it changes sign and the 
effective mass decreases monotonically from its zero temperature value with increasing temperature. 
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I. INTRODUCTION 

In the Fermi liquid theory the interacting electron sys- 
tem is composed of weakly interacting quasiparticles at 
low energies with long quasiparticle lifetimes. The ef- 
fective mass of a quasiparticle, which can be viewed as 
the bare mass of a free electron being renormalized by 
electron-electron interactions, is an important and fun- 
damental Fermi liquid parameter. For decades theorists 
have been exploring the effective mass renormalization 
in two- and three- dimensional interacting electron sys- 
tems (2DES and 3DES). In spite of this great deal of 
theoretical activity concentrating almost entirely on the 
density dependence of the effective mass renormalization, 
the temperature dependence of the effective mass has not 
been studied until very recently. Besides the consider- 
able difficulties involved in the finite temperature numer- 
ical and analytical many-body calculations in 2DES and 
3DES, the reason that this issue has not been addressed 
before can also be explained by the fact that the Fermi 
energy in 3D metals is typically lO'^K, and therefore 
any finite temperature effects are negligible. In the past 
decade, however, low density 2DES have been attracting 
attention, and several experiments have been performed 
to measure the 2D effective mass P, Q • The temperature 
dependence of the quasiparticle effective mass in 2DES is 
of considerable experimental interest since the Fermi en- 
ergy in realistic 2DES may be IK or lower, which makes 
the issue of the temperature dependence of 2D Fermi liq- 
uid parameters extremely important. In addition, the 
temperature dependence of the Fermi liquid parameters 
such as the effective mass is obviously of considerable 
fundamental theoretical significance. 

The T = quasiparticle effective mass renormalization 
in an electron system interacting through the long-range 
Coulomb interaction is one of the oldest many-body prob- 
lems in theoretical condensed matter physics, and a num- 
ber of theoretical calculations of 3D and 2D electron ef- 
fective mass have been carried out 0, 0, IE IE S H 



in the literature. In fact, the Coulomb interaction in- 
duced electron effective mass renormalization at T = is 
standard text-book material in electronic many-body 
theory. Essentially all of these calculations, both analyt- 
ical and numerical (and both 2D and 3D), are based on 
the leading-order dynamically screened interaction one- 
loop self-energy evaluation (the so-called RPA or 'GW 
self-energy approximation) because this approximation is 
really the only meaningful nontrivial calculation that can 
actually be carried out, and (perhaps more importantly) 
because this RPA self-energy is asymptotically exact in 
the weakly interacting high density regime. There have 
been a few finite temperature RPA self-energy calcula- 
tions over the years [lO, uM, U2, HEl mostly in the context 
of low dimensional systems, but none for the tempera- 
ture dependence of the effective mass renormalization in 
interacting electron systems. Very recently, Chubukov 
and Maslov considered the problem of temperature 
corrections to the 2D Fermi liquid theory for the case of a 
short-ranged interaction. In particular, they showed that 
the leading many-body temperature correction is linear 
in 2D similar to the results which we reported recently 
for the long-range Coulomb interaction in 2DES jl^ . 

In the current paper we present a calculation of the 
density and temperature dependent effective mass renor- 
malization by the Coulomb interaction in 2DES and 
3DES at arbitrary densities and temperatures (i.e. not 
necessarily restricted to high densities and low tempera- 
tures) . We work within the random-phase approximation 
(RPA), or equivalently in the ring-diagram approxima- 
tion for the self-energy, which gives exact results in the 
high-density limit (r^ <C 1) but is known to be qualita- 
tively reliable at relatively low densities as well. RPA is 
perhaps the only manageable way to perform any non- 
trivial quantitative calculations in electronic many-body 
systems, and the finite temperature RPA effective mass 
renormalization is certainly a problem of intrinsic inter- 
est. In two dimensions, our numerical results predict a 
non-monotonic effective mass temperature dependence. 
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The leading temperature dependence is linear and posi- 
tive, with the low-temperature slope being independent 
of the electron density in the high density limit. The 
temperature at which the effective mass is maximum at 
a particular density moves toward higher temperatures as 
density decreases. In three dimensions, wc find that the 
effective mass temperature dependence is non-universal 
and depends on the electron density in a complicated 
way. At very high densities, the leading correction is posi- 
tive, while at lower densities it changes sign and decreases 
monotonically from its zero temperature value. This is 
in contrast to the 2D results where the effective mass 
always increases (linearly) with temperature at low tem- 
peratures, and then decreases with temperature beyond 
a density-dependent characteristic temperature. We find 
the 3D temperature correction to the effective mass to 
be nonlinear in contrast to our 2D results. 

We express the quasiparticle effective mass m*{n, T) = 
m*{rs,T /Tp) in units of the bare band mass m (which 
is, by definition, a constant) and present our results as 
a function of the usual dimensional interaction param- 
eter rg (the average inter-electron separation measured 
in the units of Bohr radius) and the dimensionless tem- 
perature T/Tp where Tp = Ep/ks is the Fermi temper- 
ature. Note that cx n^^^^ (n^^/^) in 2D (3D) sys- 
tems, and Tp (X n {ri^/^) in 2D (3D) systems, where 
n is the appropriate 2D (per unit area) or 3D (per unit 
volume) electron density. Note that the dimensionless in- 
teraction and temperature parameters Vg and T/Tp are 
not independent parameters since they both depend on 
the electron density. We also note that Ts <C I (high- 
density) and ^ 1 (low density) limits are respectively 
the weak- and the strong-interaction limits of the electron 
system (at T = 0), and T/Tp < I and T/Tp » 1 are re- 
spectively the low-temperature (quantum) and the high- 
temperature (classical) limits. We consider the electron 
system to be a uniform jellium system with the noninter- 
acting kinetic energy dispersion being the usual parabolic 
dispersion. We use h= ks — I throughout. 

The structure of our paper is as follows: In section UTI 
we provide the formalism which we will use in this paper. 
In section Unl we explain in detail the numerical method 
we are using in the effective mass calculations. In sec- 
tion we present all our numerical results for 2D and 
3D effective mass, comparing to analytical results in the 
high density limit. In section we discuss a special ap- 
proximation method, the plasmon-pole approximation, 
and present our effective mass results using this method. 
In section IVTl we calculate the imaginary self-energy of 
quasiparticles and discuss the validity of the quasiparti- 
cle approximation at finite temperatures. We provide a 
conclusion and discussion of our results in section IVITl 



II. FORMALISM 

In this section we give the theoretical formalism, the 
basic equations, and the notations which will be used 



throughout the paper. 



A. Effective mass 

In a system of interacting fcrmions the retarded 
Green's function can be written as 



I 



eo(k) +/^ + S(k,w)^ 



(1) 



where eo(k) = fc^/2m is the spectrum of non- interacting 
fermions, fi is the chemical potential, and S(k, w) is 
the quasiparticle self-energy, the imaginary (real) part 
of which determines the lifetime (effective mass) of the 
quasiparticle. The quasiparticle energy can be obtained 
by solving the Dyson's equation 



e(k) =eo(k)+ReI](k,e(k)). 



(2) 



The quasiparticle effective mass can be written by defi- 
nition as 



m* m d 
m k dk 



k—kf 



I-|;ReS](k,^) 



I 



k dp 



^ Rcl](k,w) 



(3) 



Note that in the above equation w = is measured from 
the renormalized chemical potential /i*, which is given 
by 



H* = /_t + ReE(fcF,0). 



(4) 



All the above equations are exact, while the RPA 
approximation for S(k, w) that we are going to use is 
the first order perturbation theory in the dynamically 
screened interaction. There has been extensive discus- 
sion 0, 0, S ^1 on whether one should use exact Eq. Q 
for calculating the effective mass or it is more consistent 
to use the so-called on-shell approximation, keeping only 
the first order interaction terms in the expression for the 
effective mass (since E is calculated only to first order in 
the dynamically screened interaction): 
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™ l + fiS(k,ek)L^, 



(5) 



where = k'^ /(2m) — fi. Note that all the quantities 
on the right side of Eq. © are in the leading order in 
effective interaction. There are compelling arguments in 
favor of the latter choice: the on-shell approximation is 
believed to be more accurate as it effectively accounts 
for some higher order diagrams and satisfies the Ward 
identit y. We have extensively discussed this issue else- 
where llq. 

Obviously, the two equations for calculating m* are 
identical in the high-density limit Ts <C f. However, in 
the region of > f , they give very different results. In 
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what follows, we use Eq. ^ for all the numerical re- 
sults shown in this paper because we believe the on-shell 
approximation to be the superior one in our case. Both 
formulae give similar temperature dependence for m* (T) . 
The main qualitative results of the paper are insensitive 
to the choice of the on-shell or off-shell formula for the 
effective mass. 



B. Self-energy in the RPA approximation. 



obtained in the leading-order expansion in the dynami- 
cally screened interaction T), is conventionally called the 
RPA self-energy approximation, although the "dynam- 
ical Hatree-Fock" approximation or the "Ring Diagram 
Approximation" may be a more appropriate terminology. 

For calculations, it is more convenient to use the self- 
energy defined as a function of the real frequency co rather 
than the Matsubara one. Using the standard procedure 
of analytic continuation, one obtains the following ex- 
pression for the analytically continued self-energy: 



I]«(k,c.) 



-\-oc 



2tt 



Im G^^ (k — q, 1/ + oj) Dji (q, —v) tanh(- 



FIG. 1: Feynman diagram for RPA self-energy calculation. 
Solid lines denote the free electron Green's function and the 
dashed lines the bare Coulomb potential. 

Within RPA, the finite temperature electron self- 
energy can be expressed in terms of the Feynman dia- 
grams shown in Fig. ^ and can be written in the Mat- 
subara formalism as 7] : 

S(k, ^-n) = -T V'C/(k - q, M„ - w„OX'(q, w™), (6) 



where Vn = 7r(2n + 1)T is the fermion Matsubara fre- 
quency, ujrn = 2TTmT IS the boson Matsubara frequency 
with n and m integers, and T the temperature. The 
function I?(q, denotes the couphng to a collective 
mode (phonon, plasmon, electron-hole excitation, etc.), 
i.e., V is the bosonic propagator for the effective in- 
teraction. In our case, the function is the dynamically 
screened Coulomb interaction given by the sum of the 
ring or bubble diagrams: 



1 + 'yo(q)7r(q,cj,„)' 



(7) 



where wo(q) is the bare Coulomb interaction and 
7r(q, cjm) is the (bare) polarization operator, which is 
defined as 

7r(q,c.„) = [ ^g^'\p,.„) 

•^(")(p + q,j.„+c^„0, (8) 

where d is the dimension of the system and "(0)" denotes 
the non-interacting system. We mention that Eqs. (jZj) 
and (|HJ) together form what is called the RPA for an 
electron gas, where the bare Coulomb interaction is dy- 
namically screened by the electron dielectric function, 
which is formed from the infinite series of the polar- 
ization bubbles. The corresponding electron self-energy. 



2T 



(k — q, 1/ -I- tj) Im Dr (q, v) coth( 



2T' 



(9) 



where functions labeled with index "R" are retarded 
functions, i.e. functions analytical in the upper half- 
planes of the complex frequency. The corresponding ef- 
fective interaction can be written as: 



£'i?(q,w) 



i>(q) 



l + i;(q)nfl,(q, 



(10) 



where the retarded polarizability can be obtained from 
Eq.® using the following identities: 



nfl(q, u) = 7r(q, iujri w ir?). 



(11) 



where is a real infinitesimal positive number. 

Note that we will almost always use retarded quanti- 
ties unless otherwise stated. Thus without causing any 
confusion, we can drop the superscript "R" . 



C. Effective interaction 

The next step toward deriving the renormalization of 
mass is to obtain expressions for the effective coupling 
£)(q, Lj). We use the long-range bare Coulomb interaction 
to get 



and the effective interaction 



47re^ 

„2 ' 



Dici,Lu) = 



^ ^^o(q) 

1 + Wo(q)n(q,w) e{q,u})'' 



(12) 



(13) 



where e(q, uj) = 1 + vqH is the RPA dynamical dielectric 
function. In the RPA the full polarizability is approxi- 
mated by the bare polarizability as in Eq. ijSJl, which is 
just the bare bubble diagram. 
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Analytical properties of the propagator Hq (where TLq 
denotes the T — form for the bare polarizability with 
n denotes the finite T bare polarizability) are non-trivial 
even at zero temperature. The zero temperature polar- 
ization for 2DES and 3DES are well-known and shown 
below. For 2D T = case we have 



is defined such that 



m 
n 



2m 



3L]2 

2m' 



2m 



, (14) 



where /i is the chemical potential, the frequency oj can be 
any complex number, and the branch cut of the square 
roots are taken so that the imaginary part is positive. 
For 3D (T — 0) case we have 



nP(<z,^,M) 




where = /2ra, fi is the chemical potential and 
/i = fc^/2m = mw^/2, and the frequency uj can be any 
complex number. 

Finite temperature polarizability can be easily ob- 
tained from those at zero temperature using the following 
identity: 



oo 

n((?,w,/i;r) = J dfi' 



Uo{q,uj,fi') 
4rcosh2(4^) 



(16) 



We find Eq. to be the most convenient numerical 
method for obtaining the finite-T polarizability. 



D. Dimensionless parameters 

Our 2D and 3D electron system can be characterized 
by two parameters, namely density (n) and temperature 
(T). This immediately leads to two dimensionless pa- 
rameters Ts and T/Tp characterizing the system with Vg 
being the effective zero-temperature interaction strength 
and T/Tp being the effective temperature (note that they 
are not independent). The definition of Ts is the follow- 
ing. In 2DES rs is defined such that 



2 2 1 



kprsUB = \/2, 



(17) 
(18) 



where n is the 2D electron density, fc^ is the Fermi mo- 
mentum, and ub = (me^)"^ is the Bohr radius. In 3DES 



47rna|r3/3 = 1, 

kprsae = {9n/4)^^^ 



(19) 
(20) 



The Fermi temperature Tp = Ep — fc|,/(2m), which 
goes as Tp (X r^'^ in both 2D and 3D. 



III. NUMERICAL METHODS IN m* 
CALCULATIONS IN RPA 

In this section we explain in detail our numerical ap- 
proach for the effective mass calculation within RPA. In 
carrying out the integrations of self-energy in Eq. in 
order to obtain the effective mass, we use three differ- 
ent techniques, namely frequency sum, frequency integra- 
tion, and plasmon-pole approximation (PPA). The first 
two techniques are equivalent, and we explain them in 
detail in this section. PPA is a further approximation 
of RPAjWhich has been extensively used in the litera- 
ture [i,[l3,[l3- We discuss the PPA in section |V| Since 
there is no existing literature on the finite temperature 
effective mass or self-energy calculation to check our nu- 
merical results, it is crucial for us to use these differ- 
ent techniques to ensure the correctness of our numerical 
calculations. We mention here that our frequency sum 
results and frequency integration results agree well with 
each other. The frequency integration result is numer- 
ically relatively more noisy and therefore in this paper 
we will only show the frequency sum results. We also 
check our numerical results against the already known 
results at T = and against analytical calculations in 
the T/Tp,rs limit. 



A. Frequency integration technique 

Eq. lO gives the general formula for the RPA self- 
energy at real frequencies. It can also be written in a 
more succinct way as 



S(k,w) 



■j^vo{q)np{^^-k) 
d'^q f de 2wo('?)Inie^^(g, e) 



{2ttY J 27r e + + «ry - 
[np{£,^-\,) + nB{e)] , 



q— k 



(21) 



where np{x) — l/(exp(a;/T) + 1) is the Fermi function 
and nsix) — l/{exp{x/T) — 1) the Bose function. This 
method of calculating the self-energy involves integra- 
tion over real frequencies, and therefore we call it the 
frequency integration method. It is also known as the 
spectral or the Lehmann representation of the self-energy. 
The derivation of Eq. I|21|) from Eq. (O is given in the 
appendix. 
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The self-energy of Eq. H21() is composed of two parts: 
the exchange part and the correlation part. The (fre- 
quency independent) exchange part is also known as the 
Hartree-Fock self-energy, and its contribution to the ef- 
fective mass at T = is singular in both 2D and 3D. 
Not surprisingly, this singularity is cancelled out by con- 
tributions from the correlation part of the self-energy. 
Effective mass is derived from the self-energy through 
Eq. (0, and we therefore need to obtain the real part 
of Eq. H21|) by putting irj to be and regarding the fre- 
quency integration as a principal value integration. It is 
easy to derive from Eq. (|21|l that the imaginary part of 
the self-energy can be written as 

• K(eq-k-c^) + nF(eq-k)]. (22) 

The ImE is not needed in the effective mass calculation 
since m* is a Fermi surface property. But it is important 
to have some idea of the magnitude of ImE in order to 
ensure that quasi-particles are well defined at finite T. 

Numerically carrying out the integration in Eq. H21() 
is non-trivial: for each momentum q and frequency w, 
a three dimensional integration is required to obtain 
E(q, uj), and what makes the problem even more difficult 
is that the Ime^^(q, cj) term in the integrand is highly 
non-monotonic. A careful examination of the dynamical 
dielectric function tells us that at T = 0, Ime~-^(q, w) 
contains delta-functions at plasmon excitation frequen- 
cies, and at finite temperatures these delta-functions 
broaden into sharp peaks. Integration over these sharps 
peaks requires special care. For each q, the position (i.e., 
frequency) of the sharp peaks can be determined by solv- 
ing Re[e(q, a;)] = 0, and their weight can be determined 
from Re[e(q, 0)] using the Kramers-Kronig relations. 

One advantage of the frequency integration method is 
that in Eq. H21() we can directly put T = to obtain the 
zero temperature result, in contrast to the frequency sum 
method which we will describe in detail below. 



a transformation from the real-frequency integration to 
an imaginary frequency summation. Following the tech- 
nique of contour distortion introduced in Ref. [ij we can 
write the RPA self-energy as a sum of the Matsubara fre- 
quency along the imaginary axis: 



f dfiq 



d^q 

dS , . 

■ [?lB(Cq-k 

d'^q 



1 



1 



(27r) 



TY.v,{q) 
1 



1 



e{q,iuj„ 



1 



IUJ„ 



(^q-k - W) ' 



(23) 



where the frequency sum is over even Matsubara frequen- 
cies iujn = i2mTT with n integers. The above expres- 
sion contains three parts, namely the exchange part, the 
residue part and the line part from top to bottom in 
Eq. H23|) . The proof of the equivalence between Eq. H23|) 
and Eq. (|21|l is provided below. 

Since the exchange part exists in both Eq. (|23|l and 
Eq. 1)21(1 , we only need to consider the correlation part of 
the self-energy 



^(k,c.) 



d'^q f di^ 2vQ{q)lnie ^(g, J^) 



2t: v + UJ + irj - 



Cq-t 



(24) 



We choose the contour as in Fig. [3 It is easy to see that 
the integration over real axis can be transformed into 
integration over contour C, so that we have 



dS r dv Vo{q){t-\q,v)-l) 
[2^Y fc 2m V 



["-F(?q-k) 



- U) 

nB{v)\ . 



iV - ^q-k 

(25) 



B. Frequency sum technique 

Due to the great numerical difficulty in carrying out 
the frequency integration method introduced above (be- 
cause it involves integration over highly non-monotonic 
or singular functions), it is advisable to seek alternatives. 
At zero temperature, previous works in calculating self- 
energy and related quantities often transform the real fre- 
quency integration into integrations over imaginary fre- 
quencies using the analytic properties of the dielectric 
function. The purpose of this contour distortion is to 
avoid singularities along the real axis. At finite temper- 
ature, a similar approach can be adopted. At finite tem- 
perature, what is different from the zero temperature case 
is that we transform the integration into an imaginary 
frequency summation (or Matsubara frequency summa- 
tion) . Hu at al. [ll'l showed in detail how to perform such 



This is because e(q, w — irj) = e*(q, w -I- irj), and there- 
fore the integration of the real part of the integrand right 
above the real axis in the positive direction and right be- 
low the real axis in the negative direction cancel each 
other, and the corresponding integration of the imagi- 
nary part on these two lines are equal to each other. The 
— 1 after e~^(q, w) is inserted to make the integration 
on the arc part of contour C vanish as the radius of the 
contour approaches infinity. Now we are left to evalu- 
ate the residues within contour C, the positions of which 
are denoted by crosses in Fig. |2| Note that the analytic 
property of the dielectric function e~^(q, w) is very im- 
portant in this approach. The transformation requires 
that e~^(q, uj) is analytic in the upper and the lower half 
of the complex plane, which is true for electron gas sys- 
tems. The single residue at ^q-k ~lo — irj right below the 
real axis produced by the denominator of the integrand 
produces the residue part of the self-energy. This part 
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can be easily derived as 



over imaginary frequencies, and we have 



^res 



{2tj:Y 



1 



(26) 



The residues at a;„ = 2mrT on the imaginary axis (the 
third term in Eq. ^^), which are produced by the Bose 
function nB{v), lead to the line part of the self-energy. 
This part can be written as 



^line 



it^n - ($q-k - i^) ' 



(27) 



From Eq. ^ and Eq. ^ we have S^or = S^'C'^ + S^i^^^^ 
and we thus obtain Eq. (|23|l . 




FIG. 2: Contour of integration for the derivation of self-energy 
formula for the frequency sum method. The thick lines on real 
axis denotes the branch cut for e~^(q, tj). The crosses mark 
the poles due to the integrand; the ones on the imaginary 
axis are due to ns(aj), and the isolated pole is due to the 
denominator. 



The frequency sum method proves to be a far more effi- 
cient numerical technique for calculating the self-energy 
than the frequency integration method due to the ab- 
sence of the strong non-monotonicity and singularity in 
the real frequency dependence of the integrand. One 
thing to notice is that at high temperatures, higher 
Matsubara frequency terms can be neglected because 
(e~^(q, zwn) — 1) — > when aj„ — > oo, while at low tem- 
peratures a large number of Matsubara terms have to 
be kept in the sum in order to ensure accuracy. At zero 
temperature, the frequency sum turns into an integration 



S(k,w) 



dfiq 
df^q 



vo{q) 



va{q) 



(2^)'^e(q,^q_k-c^) 



d'^q f dv 



2ti 



{2'kY 



- 1 



(28) 



where the integration region TZi denotes the region where 
|k — q| < k, and TZ2 denotes the integration region 
where I k — q| is in between k and kp. This explicit for- 
mula for self-energy is exactly what previous works (see, 
e.g. Ref. 5j) used to calculate the zero temperature RPA 
self-energy. 

It is obvious that the frequency independent ex- 
change part of the self-energy is real. By noticing that 
e(q, — LJn) = e*(q, a;„), we can see that the line part of 
the self-energy is real as well. Thus the only contribution 
to the imaginary part of the self-energy comes from the 
residue part, which gives the same result as Eq. (|22ll in 
the frequency integration method. 



IV. RPA RESULTS FOR m'(rs,T/TF) 

In this section we present our numerical results for ef- 
fective mass in 2D and 3D electron systems within RPA. 
We first present in section llV Al results for the zero tem- 
perature effective mass to compare with earlier works. 
Our finite temperature results for 2DES are presented 
in section IIV Bl and those for 3DES in section IIV CI In 
section HV DI we present results for a model bare poten- 
tial where the Coulomb interaction is cut-off by a finite 
length so that the bare interaction is short-ranged. We 
do this in order to investigate the model dependence of 
our results. 



A. Zero temperature effective mass 

We first present our extreme low temperature results 
{T/Tp « 10-*) in Fig. El and Fig. gl to be compared 
with the existing T = results lii, i^, iSbllij iSl iSi • We 
calculate m*{rs) in the Ts = — 10 range, showing that 
the effective mass renormalization could be as large as 
4.5 for dilute {r, ~ 10) 2DES and 3 for {r, - 10) 3DES. 
We emphasize that the results presented in Fig. |3| and 
Fig. 01 are entirely based on the T — ^ limit of our finite 
temperature theory. They are in quantitative agreement 
with the existing T = 2D RPA effective mass calcula- 
tions 'Bl (which are restricted to the < 5 regime) and 
are consistent with the existing T = 3D effective mass 
calculations at low Vg . This serves as a stringent check 
on our numerical approaches. 
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FIG. 3: Calculated T ~ effective mass as a function of Vs 
in a 2DES. Inset: the result in low region. 




FIG. 4: Calculated T ~ effective mass as a function of 
in a 3DES. Inset: the result in low region. 



It is clear from Figs. 13 and 0] that both our 2D and 3D 
results show the non-monotonic dependence ofm*(rs) on 

in the high-density regime (i.e. in the <C 1 regime). 
This nonmonitonic low-rg behavior for Tn*(riAa.t T — 
has been reported in the earlier literature |3|,|5{. 

We emphasize that the numerical results given in 
Figs. 01 and H are obtained by putting T/Tp w 10""' in 
our finite-temperature formalism. 



B. Finite temperature effective mass in 2DES 

In Fig. Inland Fig. we show our calculated 2D m*{T) 
as a function of T/Tp for different values of the 2D inter- 
action parameter Ts (= 0.1 — 10). In the low temperature 
region the effective mass first rises to some maximum, 
and then decreases as temperature increases. This non- 
monotonic trend is systematic, and the value of T/Tp 
where the effective mass reaches the maximum increases 
with increasing Ts ■ The initial increase of m* (T) is almost 
linear in T /Tp as T — > 0, and the slope is almost 

independent of for very small Ts (< 1) (which is shown 
in Fig.EJ), but increases with for larger values. It is 
important to notice that this non-monotonic temperature 



dependence of m*(T) with a maximum around T /Tp < 1 
persists all the way to ^ 0, which suggests that it is 
not an artifact of our approximation scheme since RPA 
become exact as r, 0. 




FIG. 5: Calculated 2D effective mass as a function of T /Tp 
for different rg-. rs = 10 ^ 1 from top to bottom; Inset: 
Ts — b — 1 from top to bottom. Note that Tp oc rj^, making 
the absolute temperature scale lower for higher Ts values. 
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FIG. 6: Calculated 2D effective mass as a function of T /Tp 
for low Ts values: = 1.0 — + 0.1 from top to bottom. 

In Fig. [7| and Fig. |S1 we show the dependence of the 
effective mass renormalization as a function of the inter- 
action parameter Ts for a few values of fixed temperature 
(rather than fixed T/Tp, remembering that Tp cx rj^ 
since Tp cx n and Vs oc n""'^/^). Fig. [7|shows the effective 
mass for high T and large values while Fig. |S| concen- 
trates on the low T region. The calculated m*{rs) for 
fixed T values are quite striking: For low fixed values 
of T, m* /m initially increases with even faster than 
the corresponding T = result, eventually decreasing 
with Ts at large enough values (where the corresponding 
T/Tp values become large enough). This nonmonotonic 
behavior of m*(rs) as a function of for fixed tempera- 
tures showing a temperature-dependent maximum (with 
the value of at which the m* peak occurs decreas- 
ing with increasing T as in Fig. [T)) is complementary to 
the nonmonotonicity of m*{T) in Fig. [S] as a function of 
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T/Tp (at fixed Ts) and arises from the relationship be- 
tween the dimensionless variables T/Tp (oc r~'^) and 

— 1/2 

(cx Tp ) due to their dependence on the carrier den- 
sity (i.e. Tp oc n and oc n~^/'^). At large and 
high temperature, Fig. shows that the effective mass 
increases from below unity with increasing r^. This is 
the region where the exchange part of the self-energy 
dominates, and it can be easily shown that the exchange 
self-energy produces this peculiar effect on the depen- 
dence of m*(rs) at fixed high T values. Since this region 
is hardly accessible by experiments, and moreover the 
quasiparticles may not even be well-defined at such high 
T /Tp values, we do not further discuss the physics re- 
lated to this region. 



consistent with recent experimental findings 1] , our the- 
oretical calculations predict (according to Figs. |S1 and [7| 
as modified by subband form factors) m*/m to be en- 
hanced by a factor of 2 — 4 for the experimental densities 
and temperatures used in recent measurements 1]. Due 
to the approximate (i.e. RPA) nature of our theory we 
do not further pursue the comparison with experimental 
data in this paper since the main goal of this paper is 
to discuss the temperature dependence of m*{rs,T/Tp) 
which has not yet been reported in the literature. A di- 
rect experimental observation of an increasing m* (T) at 
low temperatures in 2DES will be a striking confirmation 
of our theory. 




FIG. 7: Calculated 2D effective mass as a function of at 
fixed value of temperatures. T is in the unit of Tp at = 1. 




FIG. 8: Calculated 2D effective mass as a function of at 
fixed value of temperatures. T is in the unit of Tp at = 1. 
This plot is similar to Fig. but concentrating on the low 
temperature region. 

One immediate consequence of our results shown in 
Figs. El and His that m*{T/Tp,rs) = m*(T,n) in 2DES 
could show a strong enhancement at low (but finite) tem- 
peratures and low electron densities (large r^). Compar- 
ing with the actual system parameters for 2D electrons 
in Si inversion layers and GaAs heterostructures (and 
taking into account the quasi-2D form factor effects 
neglected in our strictly 2D calculation) we find that, 



C. Finite temperature effective mass in 3DES 

In Fig. inland Fig. ^ we show our calculated 3D m*{T) 
as a function of T/Tp for different values. In Fig. 
Ts varies from 1 to 10 while in Fig. 1101 is from 0.1 
to 1. The 3D temperature dependence of the effective 
mass shows very different characteristics from that of 2D. 
Fig. [51 shows that for > 1 the effective mass decreases 
monotonically with increasing T at low temperatures. 
However for Vg « 1, as shown in Fig. ^| the effec- 
tive mass increases with increasing T in the temperature 
region we are considering. We therefore conclude that in 



3DES the sign of the slope 



d{m* / m) 



is non-universal, 



d{T/TF)\T = 

which differs from that of 2DES where the above men- 
tioned slope is always positive for all r^. Another in- 
teresting feature is that the sign of '^j^^/^j \t=o matches 
the sign of — '^^^"^^"'^ \t=o very well. In particular, m*{T) 
decreases with increasing T (at low T) in the "larger" 
r^-regime where the corresponding T = m*{rs) shows 
an increasing mass with increasing r^. Similarly, m*{T) 
increases (at low T) with increasing T in the rg-regime 
where the corresponding m* (r^ ; T = 0) shows decreasing 
m* with increasing r^. 
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FIG. 9: Calculated 3D effective mass as a function of T/Tp 
for different rg-. rs = 10 ^ 1 from top to bottom. 
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FIG. 10; Calculated 3D effective mass as a function of T/Tf 
for low fs values: Ts = 1.0 — > 0.1. 



D. Model short-range bare interaction 

So far in all of our calculations we have been using the 
realistic long-ranged Coulomb interaction for the bare 
potential as in Eq. (|12|l . A question naturally arises: 
how is the temperature dependence that we find in our 
calculations related to the long range nature of the in- 
teraction between electrons? Therefore we also calculate 
the effective mass in 2DES and 3DES using simple (pa- 
rameterized) finite-range interaction model 



q + a 

+ a? 



(29) 



where a is the cut-off wavevector which eliminates the 
long wavelength Coulomb divergence. 

Our numerical calculation shows that as a/kp 0, we 
recover the m* (T) behavior of the bare Coulomb inter- 
action results in both 2D and 3D. As a/kp increases, the 
mass renormalization in both 2D and 3D is suppressed, 
but all the qualitative features of the temperature de- 
pendence persist. In 2DES, as a/kp increases, the tem- 
perature where the effective mass reaches the maximum 
decreases, and the effective mass enhancement (from the 
T = value to the maximum) decreases, but the linear- 
T dependence at low T and the non-monotonic trend 
remain unchanged. In 3DES, as a/kp increases, the Vs 
region where \t=o < shrinks, but the consis- 

tency between the sign of — '^^'^(^J^^ \ t=o and the sign of 

d{T/TF) \T=o remams. 

From these result we conclude that the qualitative fea- 
tures of the temperature dependence are model inde- 
pendent and not peculiar to the bare interaction being 
Coulombic. This conclusion is further reinforced by the 
recent report of a linearly T dependent electronic specific 
heat in a short-range interaction model \T^ . It may be 
worth while, however, to note that RPA is specific to the 
long-range Coulomb interaction in giving an exact result 



in the high-density — > limit, and there is nothing 
special about RPA in the case of short-range interaction. 



V. PLASMON-POLE APPROXIMATION 

We now apply a simple-to-use dynamical approxima- 
tion to calculate m*(T). The plasmon-pole approxima- 
tion (PPA) has often been used [a, 113 to obtain the 
electron self-energy in the literature. It is a simple tech- 
nique for carrying out the frequency sum or integration in 
the RPA self-energy calculation by using a spectral pole 
(i.e. a delta function) ansatz for the dynamical dielectric 
function e(k, w). In other words, it is an approximation 
to the RPA. The PPA ansatz assumes that 



-2Im 



1 



e(k,w) 



Ck{Siio-Ldk)^6iLu + Qk)), (30) 



where the pole Cjk and the spectral weight Ck of the PPA 
propagator in Eq. H30() are determined by using the the 
Kramers-Kronig relation (i.e. causality) 



Re 



1 



e{k,0) 



2 

1 + - 

TT 



1 , 1 

u) e{k,uj) 



and the /-sum rule (i.e. current conservation) 

1 



uidujlm- 
^[k,uj) 



ulik). 



Putting Eq. ^ in Eq. (EH and Eq. ^ we have 
Ck - 7rcjp(fc)v/l-Ree-i(fc,0), 

LUp{k) 



^1 -Ree~i(fc,0)' 



(31) 

(32) 

(33) 
(34) 



where ujp{k) in Eqs. and is the long wave- 

length plasmon frequency which is defined as 



lim Re[e{k,uj)] = 1 - 

It is well-known that in 2DES 

27rne^ . 



ul{k) 



uj'pik) 



and in 3DES 



,(fc) = 



Aime 



(35) 



(36) 



(37) 



We mention that Qk in Eq. 1)30(1 does not correspond 
to the real plasmon dispersion in the electron liquid, but 
simulates the whole excitation spectra of the system be- 
having as an effective plasmon at low momentum and as 
the single-particle electron-hole excitation at large mo- 
mentum, as constrained by the Kramers-Kronig relation 
and the /-sum rule. Details on the PPA are available in 



10 



literature @J0,^3|, including its finite-temperature gen- 
eralization [l2|- The PPA, which is known to give results 
close to the full RPA calculation of self-energy, allows a 
trivial carrying out of the frequency sum in the retarded 
self-energy function leading to: 



Rel](k,w) 



_fq_ 
+ 



Wo('7)"-F(Cq-k) 



(e 



q-k 



(38) 



where Cq and Uk only depend on e(fc, 0) at finite tem- 
peratures, and are determined by Eq. (|33|) and Eq. (|34ll . 
Obviously the PPA provides a great simplification of the 
problem since the most numerically demanding part of 
the calculation (the frequency sum or integration) is triv- 
ially done. It should be noted, however, that although 
the PPA is known to produce a reliable approximation 
to ReE, it, by definition, fails completely for ImS. 

We present our PPA results for the 2D effective mass 
as a function of T/Tp at fixed values in Fig. ^] One 
immediate observation by comparing Fig. [S] and Fig. 1111 
is that even though PPA provides a very good approx- 
imation for the self-energy (indeed, our numerical re- 
sults for PPA self-energy and RPA self-energy match very 
well), it fails to provide accurate result for the effective 
mass. The zero temperature effective mass generated by 
PPA is almost half of that from RPA, and the tempera- 
tures where m* maximizes shift to higher T values in the 
PPA compared with RPA. But the qualitative behavior 
of m*{rs,T/TF) is similar in the PPA and RPA for the 
2DES as is clear by comparing Figs. 1111 and 1^ 



r = 10 




FIG. 11: Calculated 2D PPA effective mass as a function of 
T/Tf at fixed value of 7-3. 

From our results of 3D PPA effective mass calculation 
presented in Fig. El we can see that they are different 
from RPA results even qualitatively. In fact, our RPA 
results for m*{rs,T/Tp) are similar in both 2D and 3D. 
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FIG. 12: Calculated 3D PPA eflective mass as a function of 
T/Tf at fixed value of Vs- 



VI. QUASIPARTICLE DECAY 

The quasiparticle decay rate (or the inverse lifetime) 
is given [l^ by the imaginary part of the self-energy. 
As we have discussed in section BlI Bl the imaginary part 
of the quasiparticle self-energy can be calculated from 
Eq. I|23() . It is also obvious that only the second term 
in Eq. (|23|l contributes to the imaginary self-energy: the 
first term is obviously real, and the last term is also real 
because e(q, — Wn,) = e*(q, ciJn). Thus we have 



ImE(k,w) 



-—^vo[q)lrn—— 

i^-^r e(9,Cq-k) 

"-B(Cq-k ?T-F(Cq-k)] • (39) 



Fig. El and Fig. El show the calculated imaginary 
self-energy on the Fermi surface in 2D and 3D respec- 
tively. The quasiparticle decay (i.e. finite ImE) here 
arises entirely from having a finite temperature. The 
results show that the magnitude of the imaginary self 
energy increases with increasing rg and T/Tp. It is ob- 
vious from Eq. (|39|l that the imaginary self-energy van- 
ishes on the Fermi surface at T = as it must since 
the quasiparticles are perfectly well-defined at T = for 
k = kp- As T increases, the magnitude of imaginary 
self-energy remains small compared to the Fermi energy 
up to a certain temperature, and the quasiparticles on 
the Fermi surface remain well-defined up to that temper- 
ature. The important question is whether the finite tem- 
perature quasiparticles are sufficiently well-defined for 
the interesting behavior of m*(T) we discussed in sec- 
tion llVl to be experimentally observable. If the quasipar- 
ticles are ill-defined (i.e. ImE(fci?) > Ep in the tempera- 
ture regime of interest) then obviously all the interesting 
temperature dependence of m*(T) predicted by us is only 
of academic interest since the large broadening will make 
it impossible to define quasiparticles, let alone their ef- 
fective mass. Examining the results of Figs. 1131 and 1141 
compared with those presented in section IIVI it is clear 
that there is a well-defined regime of {rs,T/Tp) values 
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where m* (T /Tp) shows non-trivial temperature depen- 
dence with the condition Ep ^ |IniE(fci?)| weh-satisfied 
so that quasiparticles are weh-defined. Although this 
is not unexpected since |IniE(r)| - for T/Tp < 1 
whereas m*{T)/m 1 + 0{T) in 2D, it is nevertheless 
important to see that ImE remains small in magnitude 
in the (rs,T/Tp) regime of interest. 




FIG. 13: Calculated magnitude of the 2D RPA imaginary 
self-energy of quasiparticles on Fermi surface as a function of 
Ts at different values of T/Tp- 





9 




0.4 








0.8 




8 




1.2 




7 




1.6 

2.0 








2.4 




6 




2.8 








3.2 




5 




3.6 








4.0 




4 




w 








3 






FIG. 14: Calculated 3D RPA imaginary self-energy of quasi- 
particles on Fermi surface as a function of at different values 
of T/Tf. 

Earlier theoretical work on the quasiparticle damp- 
ing of 2D interacting electron systems can be found in 
Refs. HllISIll. 



VII. DISCUSSION AND CONCLUSION 

In this work we have obtained detailed results for 
the temperature dependence of the quasiparticle effective 
mass, m*{rs,T/Tp), at arbitrary values of temperature 
and density in 2D and 3D electron systems interacting via 
the long range Coulomb interaction. Our central approxi- 
mation is the RPA (i.e. the dynamically screened Hatree- 
Fock self-energy approximation), which is the leading- 
order one-loop self-energy calculation in a dynamically 



screened effective interaction expansion. RPA is exact in 
the high-density (r^ 0) limit at T = 0, and is there- 
fore a controlled nontrivial approximation which is em- 
pirically known to work well for > 1 (e.g. metals with 

~ 3 — 6 and 2D semiconductors with 1 — 10). 

We also calculate the finite-temperature imaginary self- 
energy (i.e. the quasiparticle decay rate or broaden- 
ing) to ensure that the broadening remains small in the 
{rs,T/Tp) parameter regime of our interest where m*{T) 
shows interesting temperature dependence. 

As mentioned earlier in the paper, it is well-known that 
at T = 0, m*{rs) can be exactly calculated (in both 2D 
and 3D) in the asymptotic ^ limit by systematically 
expanding the RPA self-energy since ring diagrams (in- 
cluded in the RPA) are the most divergent diagrams in 
the — > limit. Such a zero temperature Tg-expansion 
of RPA gives the following formula for m* (r^) in both 2D 
and 3D: 



m*{rs) 



l + ars{b + \nrs)+0{rl), (40) 



where a, h are constants of order unity. What we find 
numerically is that the leading temperature correction 
to this effective mass formula is linear in T/Tp in 2D 
and nonlinear in T /Tp in 3D. In this paper we have cal- 
culated m* (tsiT /Tp) numerically for the one-loop dy- 
namically screened Hatree-Fock RPA self-energy theory 
for arbitrary and T /Tp finding nontrivial temperature 
dependence of the effective mass at all densities. 

Our most important result is the unexpected discov- 
ery of a strong temperature-dependent quasiparticle ef- 
fective mass m*(T) at low temperatures in 2DES. Since 
the temperature scale for the temperature dependence of 
m*(T) is the Fermi temperature which tends to be high 
(~ 10"* if) in the 3D electron liquids (i.e. metals), our 
temperature-dependent effective mass results for 3D sys- 
tems are mostly of theoretical interest since any actual 
T-dependence of m*{T) in the T/Tp < lO""^ regime will 
be miniscule. Our numerical results for the calculated 
m*{T) in 2D systems are consistent with a Zmear leading- 
order temperature correction for the 2D quasiparticle ef- 
fective mass: Results in Figs.|51and|Slcan be well fitted to 
the formula m*(T) « 1 + A^^ {r^) + B"^^ {rs){T /Tp) + ■ ■ ■ 
for small T/Tp where the slope B^^{rs) seems to be a 
constant independent of (i.e. density) at least in the 
high-density (r^ ^ 1) limit; for > 1 the slope B^^{rs) 
has a weak density dependence increasing somewhat with 
increasing (but our approximation scheme, RPA, be- 
comes less quantitatively reliable at large r^, therefore it 
is possible that the slope d{m* / m) / d{T /Tp) is indeed in- 
dependent of Ts in the T ^ limit). In addition to this 
interesting (and unexpected) linear leading-order tem- 
perature correction to the quasiparticle effective mass, 
we also find B^^{> 0) to be positive for all rg, indicting 
that in 2DES, the leading-order temperature correction 
to the effective mass is positive. Thus, m*{T) increases 
with increasing T at first before eventually decreasing as 
T/Tp increases substantially, leading to a maximum in 
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m*{T) at some intermediate temperature T*{rs) ^ 0.5Tp 
which is only weakly density dependent (except of course 
through Tf itself). All three of these 2D findings (i.e. lin- 
ear leading-order T/Tp dependence of m*, increasing m* 
with T/Tp at low temperatures, and the nonmonotonic 
behavior with a maximum in m*{T/Tp) occurring at 
T* ~ 0.5Tp) are surprising and unexpected. In principle, 
these predictions can be experimentally tested since our 
calculations presented in section FVll show that the quasi- 
particles remain reasonably well-defined (i.e. the broad- 
ening, ImS, remains small) all the way to T* and perhaps 
even above T* . This is reasonable since the many-body 
correction to m* is linear in T/Tp whereas the broad- 
ening ImE ~ [T/Tpf, ensuring that for T/Tp < 1, the 
quasiparticle effective mass is a well-defined quantity. In 
contrast to the linear (with positive slope) leading-order 
T-dependence we find for all in our calculated 2D 
m*{T), our 3D results show non-universal m* {rs,T /Tp) 
behavior. In 3D, m*{T/Tp) increases with increasing 
T/Tp at low temperatures only for very high densi- 
ties (small rs) - for larger Ts values m*{T) decreases 
monotonically with increasing temperature (in sharp con- 
trast to the striking non-monotonicity in m*{T) in 2D) 
and this decrease is more consistent with a nonlinear 
leading-order temperature dependence (rather than a lin- 
ear one as in 2D). Our best guess for our numerical re- 
sults shown in Fig. |51 and ^| is the following equation: 
m*/m K l-VA^^{rs)+C^^{T/Tpy \n{Tp/T) + - ■ ■ , where 
I is a number of the order one (note that numerically fix- 
ing the number I is difficult and need much more work), 
> for Ts < r* and < for > r* with r* be- 
ing approximately the Vg value where A^^ changes from 
being negative to positive. 

We comment that our numerical results for to* (T) are 
consistent (actuall y a g ree very well) with the very recent 
analytical work on the temperature corrections to 

the effective mass renormalization in 2D and 3D Fermi 
liquid. The analytical work is necessarily restricted to 
the Ts — > and T/Tp — s- limit where the infinite series 
of ring diagrams to the electron self-energy (depicted in 
Fig-ffl) provides an exact leading-order asymptotic answer 
to the problem with the following result 

^!:(r^.i + .4(,, + B,..,f^) 

TO \Tp J 

with B{rs) = B^D, a constant, in 2D, and B{rs) = 
in 3D. Our numerical results are consistent with this ex- 
act result, but our numerical results apply also in the 
non-asymptotic region where T/Tp and are not neces- 
sarily small. In this non-asymptotic regime (where is 
not small, actually may be large in 2D semiconductor 
systems) RPA is by no means an exact theory, but we 
have recently argued that RPA remains extremely 
well- valid (if somewhat uncontrolled) even for ^ 1. 
We also emphasize a point in this context that seems 



not to have been widely appreciated in the literature. 
The point is that RPA becomes a progressively better 
approximation as T/Tp increases at a fixed (for any 
rs), because the system is becoming more classical in the 
T/Tp ^ 1 regime where RPA is again exact. Thus in the 
{rs,T/Tp) parameter space (see Fig. llSf) RPA is exact as 
— > (the high density limit) and as Tp — > (the high- 
temperature or equivalently the low density limit) or as 
T oo. Thus the regime of validity of RPA is greatly 
enhanced at finite temperature, and in fact even at very 
large (i.e. very low density) RPA becomes exact as T 
is raised (because T/Tp ^ 1 limit is more easily achieved 
at low densities). 
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FIG. 15: Schematic validity of RPA. The shaded area denotes 
the region where RPA is considered to be valid. Line A de- 
notes a certain density value above which RPA is valid at 
T = (e.g. the vertical line A may correspond to the rs — 1 
condition so that for higher density, i.e. to the right to line 
A, RPA is valid even at T = 0). Line B denotes the line 
of Tp oc n. In the region above line B, RPA is again valid. 
Therefore for any fixed value of density n (or equivalently 
fixed rs), RPA is a better approximation with increasing T, 
whereas for fixed value of temperature, low density values (or 
large r,) counter- intuitively makes RPA to valid again since 
RPA is valid for T >Tf. 

Finally, we comment on the anomalous (often referred 
to as "nonanalytic" [l4|) nature of the temperature cor- 
rections to the quasiparticle effective mass in 2D systems 
(but not in 3D) as manifested in the linear leading-order 
temperature correction we find in interacting 2D electron 
systems. This particular feature is apparently generic 
in 2D and not due merely to our using the long-range 
bare Coulomb interaction, because in Ref. 14] the same 
linear-T correction is found in calculations using a zero- 
range bare interaction although the sign of the slope is 
negative in the zero-range interaction case. This kind 
of (leading-order) linear temperature correction is quite 
common in 2D electron systems due to the peculiar form 
of the 2D polarizability with a T = cusp at 2kp. This 
leading-order linear-T correction is interesting because 
the naive expectation in a Fermi system (based on the 
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usual Sommerfeld expansion of the Fermi functions) is 
that the leading-order correction in a "normal" situation 
should always be 0{T/TfY for electronic properties. 
In 2D electron systems it seems that the generic situation 
is "anomalous" , i.e. the leading-order temperature cor- 
rection is 0{T/Tp) rather than the "normal" quadratic 
Fermi behavior expected on the basis of analytic Som- 
merfeld expansion of Fermi functions. In 2D interacting 
electron systems, therefore, all leading-order thermal cor- 
rections to electronic properties are much stronger (by 
a factor of Tp/T, which is a large number as T ^ 0) 
than the quadratic Fermi gas behavior. This anomalous 
nonanalyticity, which may have important consequences 
for fermionic quantum critical phenomena, obviously has 
important experimental implications since it is much eas- 
ier to observe a linear temperature correction than a 
quadratic one at low temperatures. 

This work is supported by the US-ONR, the NSF, and 
the LPS. 



Using Kramers-Kronig relations for e ^(q, v) in the above 
equation, we have 
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Appendix 

Here we provide a proof of for equivalence between 
Eq. 123 and Eq. PT|). 
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